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We construct explicit examples of globally regular static, spherically symmetric solutions in general 
relativity with scalar and electromagnetic fields which describe traversable wormholes (with flat 

. and AdS asymptotics) and regular black holes, in particular, black universes. A black universe is a 

nonsingular black hole where, beyond the horizon, instead of a singularity, there is an expanding, 

■ asymptotically isotropic universe. The scalar field in these solutions is phantom (i.e., its kinetic 

energy is negative), minimally coupled to gravity and has a nonzero self-interaction potential. The 
configurations obtained are quite diverse and contain different numbers of Killing horizons, from 
zero to four. This substantially widened the list of known structures of regular BH configurations. 
Such models can be of interest both as descriptions of local objects (black holes and wormholes) 
' and as a basis for building nonsingular cosmological scenarios. 
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1 Introduction 



One of the basic problems of black hole (BH) physics is the existence of curvature singularities 
beyond the event horizons in the well-known Schwarzschild, Reissner-Nordstrom, Kerr and other 
solutions of general relativity and their analogs in other metric theories of gravity. For full under- 
standing of BH physics and geometry it is highly desirable to get rid of singularities, and this is 
usually connected with the hopes for a future quantum gravity. Still of great interest are attempts 
to construct non-singular BH models in the framework of classical gravity, and different classes of 
such models have been described in the literature. One of such classes, termed black universes [1,2], 
is in our view of particular interest since it combines the properties of wormholes (no center, and a 
regular minimum of the area of coordinate spheres), BHs (a Killing horizon separating static and 
non-static space-time regions) and non-singular cosmological models (at large times the non-static 
region reaches a de Sitter mode of isotropic expansion). The black universe models make possible a 
cosmological scenario where a phantom-dominated gravitational collapse in some "mother" universe 
creates our Universe whose expansion begins from a horizon, and the next stages are isotropization 
and de Sitter inflationary expansion. Other kinds of regular BHs discussed in the literature are 
classified in [2]; see also the conclusion of the present paper. 

In the models described in [1,2], the material source of gravity is a phantom scalar field that 
differs from the canonical one by the sign of its kinetic energy. Such a field has been repeatedly 
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discussed as a possible supporter of wormhole geometries ( [3,4] and a great number of later papers; 
see [5,6] for recent reviews). The possible existence of phantom fields in the Nature is to a large 
extent favored by modern cosmological observations, indicating that the accelerated expansion of our 
Universe may be caused by a dominating dark energy density with the pressure to density ratio w 
smaller than -1. One can note that values w < — 1 seem to be not only admissible but even preferable 
for describing an increasing acceleration, as follows from the most recent estimates: w = — 1.10±0.14 
(la) [7] (according to the 7-year WMAP data) and w = -1.0691q:o92 [ 8 ] (mainly from data on 
type la supernovae from the SNLS3 sample). In this connection, cosmological models with phantom 
scalar fields, i.e., those with a negative kinetic term, have gained considerable attention in the recent 
years (see, e.g., [9,10] and references therein). There are theoretical arguments both pro et contra 
phantom fields, and the latter seem somewhat stronger, see, for instance, a discussion in [11]. 

In this paper, accepting the existence of a phantom scalar as a working hypothesis, we would like 
to discuss new features of wormhole and black-universe configurations which appear if, in addition 
to a scalar field, an electromagnetic field is invoked as a source of gravity. As in [1,2], we deal with 
static, spherically symmetric space-times, therefore the only kinds of electromagnetic fields are a 
radial electric (Coulomb) field and a radial magnetic (monopole) field. It should be stressed that in 
the latter case it is not necessary to assume the existence of magnetic charges (monopoles): in both 
wormholes and black universes a monopole magnetic field can exist without sources due to a specific 
space-time geometry. In the wormhole case it perfectly conforms to Wheeler's idea of a "charge 
without charge" [13]: electric or magnetic lines of force simply thread the wormhole. In the case of 
a black universe, the picture is different on different sides: in the static region a possible observer 
sees a black hole with an electric or magnetic charge; in the cosmological region, this corresponds 
to a homogeneous primordial electric or magnetic field. For definiteness, we will speak of magnetic 
fields. 

One of motivations for the present study was that modern observations testify to a possible 
existence of a global magnetic field up to 10~ 15 Gauss, causing correlated orientations of sources 
remote from each other [12], and some authors point out the possible primordial nature of such a 
magnetic field. 

The paper is organized as follows. In Section 2 we present the basic equations and make some 
general observations. In Section 3 we obtain explicit examples of wormhole, black-universe and other 
regular black hole solutions using the inverse problem method. Section 4 contains a discussion and, 
in particular, some numerical estimates concerning the possible magnetic field strength at different 
stages of the cosmological evolution. 



2 Basic equations 

We consider the action 4 

r^-gofx \R + 2eg^d^d u <P - 2V{<j>) - F^F^ 
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(1) 



where R is the scalar curvature, g = det((? At!/ ), and F^ is the electromagnetic field tensor, e — +1 
corresponds to a normal scalar field and e — — 1 to a phantom one. 

The general static, spherically symmetric metric can be written in the form 

ds 2 = A(u)dt 2 -4^~ r 2 (u)dil 2 , (2) 



A(u) 



4 We choose the metric signature (- 
the energy density. 



) , the units c = h = 8ttG = 1 , and the sign of such that T§ is 
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where we are using the so-called quasiglobal gauge googu = — 1; A(u) is called the redshift function 
and r(u) the area function; dVl 2 = (d6 2 + sin 2 9 dip 2 ) is the linear element on a unit sphere. The 
metric is only formally static: it is really static if A > , but it describes a Kantowski-Sachs type 
cosmology if A < 0, and u is then a temporal coordinate. In cases where A changes its sign, regions 
where A > and A < are called R- and T-regions, respectively. 

Let us specify which kinds of functions r(u) and A(u) are required for the metric (2) to describe 
a wormhole or a black universe. 

1. The range of u should be u G M., where both A(u) and r(u) should be regular, r > 
everywhere, and r — > oo at both ends. 

2. A flat, de Sitter or AdS asymptotic behavior as u — >■ ±oo. 

3. In the wormhole case, absence of horizons (zeros of A(u)), and flat or AdS asymptotics at 
both ends. 

4. In the black-universe case, a flat or AdS asymptotic at one end and a de Sitter asymptotic at 
the other. 

The existence of two asymptotic regions with r ~ \u\ (by item 2) requires at least one regular 
minimum of r(u) at some u — uo, at which 

r = r > 0, r' = 0, r" > 0, (3) 

where the prime stands for d/du. (In special cases where r" = at the minimum, we inevitably 
have r" > in its neighborhood.) 

The necessity of violating the weak and null energy conditions at such minima follows from the 
Einstein equations. Indeed, one of them reads 

2Ar"/r = -(T t t -T:), (4) 

where T£ are components of the total stress-energy tensor (SET). 

In an R- region (A > 0), the condition r" > implies T\ — < 0; in the usual notations T\ = p 
(density) and — = p r (radial pressure) it is rewritten as p + p r < , which manifests violation 
of the weak and null energy conditions. It is the simplest proof of this well-known violation near a 
throat of a static, spherically symmetric wormhole ( [14]; see also [6]). 

However, a minimum of r(u) can occur in a T- region, and it is then not a throat but a bounce 
in the evolution of one of the Kantowski-Sachs scale factors (the other scale factor is [— A(-u)] 1 / 2 ). 
Since in a T-region t is a spatial coordinate and u temporal, the meaning of the SET components 
is — T\ = p t (pressure in the t direction) and = p; nevertheless, the condition r" > applied 
to (4) again leads to p + pt < , violating the energy conditions. In the intermediate case where 
a minimum of r(u) coincides with a horizon (A = 0), the condition r" > holds in its vicinity, 
along with all its consequences. Thus the energy conditions are violated near a minimum of r in 
all cases. 

In what follows, we will assume that the space-time is asymptotically flat as u — > oo and consider 
different behaviors of the metric as u — > — oo. 

The scalar field <f>(u) involved in the action (1) in a space-time with the metric (2) has the SET 

T;[s] = eAiuWiu) 2 diag(l, -1, 1,1)+ 5;V(u). (5) 

The solutions of interest for us correspond to e = — 1 but we preserve both values of e in the 
equations for generality. 
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The electromagnetic field compatible with the metric (2) can have the following nonzero com- 
ponents: 

F i = — Fio (electric), and F 23 = -F 32 (magnetic), 
such that 

F 01 F 01 = -q 2 Jr\u), F 23 F 23 = q 2 Jr\u), (6) 

where the constants q e and q m have the meaning of electric and magnetic charges, respectively. 
The corresponding SET is 

7£[e] = -^diag(l, 1, -1, -1), q 2 = q 2 e +ql (7) 
r {U) 

Thus the electromagnetic field equations have already been solved in a general form, and we are 
left with the set of Einstein and scalar field equations. It can be written as follows: 

2{Ar 2 <f>')' = er 2 dV/d<f), (8) 

(A'r 2 )' = -2r 2 V + 2q 2 /r 2 , (9) 

r"/r = -e<f)' 2 , (10) 

A{r 2 )" - r 2 A" = 2 - 4q 2 /r 2 , (11) 

- 1 + A'rr' + Ar' 2 = r 2 (eA0' 2 - V) - q 2 /r 2 , (12) 

The scalar field equation (8) follows from (9)-(ll), which, given the potential V ((/>), form a deter- 
mined set of equations for the unknowns r(u), A(u), 4>(u) . Eq. (12) (the (j) component of the 
Einstein equations), free from second-order derivatives, is a first integral of (8)— (11) and can be 
obtained from (9)— (11) by excluding second-order derivatives. Moreover, Eq. (11) can be integrated 
giving 

^'(")-'' 4 (^)' = -2- + ^/-^ y . (13) 
where B{u) = A/r 2 . 

Let us note that Eqs. (8)-(12) in the case of a massless scalar field <p have been solved long ago, 
in [15] for e = +1 and in [3] for e — — 1. At e = +1 all such solutions possess a central singularity; 
with a phantom scalar (e = —1), there are both singular solutions and twice asymptotically flat 
wormholes [3] but nothing like black universes. 

We here seek solutions with a nonzero potential V ((/>). It is known [16] that Eqs. (8)-(12) lead 
to a very narrow choice of possible global space-time structures in the case q — 0. Indeed, due to 
(11), if q = 0, the function B{u) cannot have a regular minimum, therefore it can have at most 
two zeros (which coincide with zeros of A(u) and hence correspond to horizons), and if the model 
is asymptotically flat, say, at large u, only a single simple horizon is possible. We shall see how a 
nonzero charge q changes the situation. 
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3 Some particular models 
3.1 Solutions 

If one specifies the potential V(4>), it is, in general, very hard to solve the field equations. Alter- 
natively, to find examples of solutions possessing particular properties, one may employ the inverse 
problem method, choosing some of the functions r(u), A{u) or <f)(u) and then reconstructing the 
form of V{4>) . We will do so, choosing a function r(u) that can provide wormhole and black-universe 
solutions. Given r(u) and the charge q, the function A(u) is found from (13) and V(u) from (9). 
Furthermore, 4>(u) is found from (10), and, as long as r"/r ^ 0, we obtain a monotonic function 
4>{u) which then yields an unambiguous function V{4>). 

A simple example of the function r(u) compatible with the requirements 1-4 is [1] 



r(u) = Vu 2 + b 2 = bVx 2 + 1 (14) 

where x = u/b, and b > is an arbitrary constant (the length scale). Evidently, r"(x) > 0, as 
required, and thus we automatically fix e — —1; we also have r b\x\ at large \x\ . 

Let us formally put b = 1 , which will actually mean that the length scale is arbitrary but the 
quantities r, q, m (the Schwarzschild mass in our geometrized units) etc., with the dimension of 
length, are expressed in units of b, the quantities B, V and others with the dimension (length) -2 
in units of b~ 2 , etc.; the quantities A and are dimensionless. 

Now, the expression for B' = B'(x) can be written as 

2 

B '( x ) = 71~, — 5vj(P ~ x + 2g 2 arctanx), (15) 
(1 + x z ) z 

where p is an integration constant. Further integration gives 

^/ n ^ l + q 2 +px ( 2q 2 x \ 9 9 

B(x) = B H — ^— + I p H — - 1 arctanx + q 2 arctan 2 x, (16) 

1 I tjC V 1 I tic J 

where Bq is one more integration constant. 

Now suppose that our system is asymptotically flat at x — >■ +oo. Since B = A/r 2 and A — > 1 
at infinity, we require -B — >■ as rr — > oo and thus fix B as 

5 = -irp/2 - n 2 q 2 /A. (17) 

Furthermore, comparing the asymptotic expression A = 1 — 2m jx + o(x) for A(x) with what is 
obtained from our expression for A = Br 2 , we find a relation between the Schwarzschild mass m 
and our parameters p and q: 

p = 3m — irq 2 . (18) 

Thus B is a function of x and two parameters, the mass m and the charge q . 

Now we know the metric completely, while the remaining quantities <f>(x) and V((f>(x)) are easily 
found from Eqs. (10) and (9), respectively: 

0(x) = ± arctanx + 0o, 0o = const; (19) 

q 2 1 r 

V(x) = r-7T ; 2q 2 (3x 2 + 1) arctan 2 x 

v ; (1 + x 2 ) 2 2(1 +x 2 ) L y v ; 

+ (18x m — 67rx g + 12xg — 2rcq + 6m) arctan re 

+ g 2 (^V - 67ra; + ^vr 2 + 6^ - m^nx 2 - 18x + 3tt)] . (20) 
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Figure 1: Plots of B(x) (left) and V[x) (right) for symmetric configurations. Curves Al, A2, A3, A4 
correspond to q = 0.7, 0.77, 0.82461, 0.9, respectively. The lower panel shows the Carter-Penrose diagrams 
of the corresponding space-times. 

Thus has a finite range: G (— 7r/2,7r/2) (putting O = without loss of generality), which 
is common to kink configurations. We also have x = u/b = tan0, whose substitution into the 
expression for V(x) gives V((j)) defined in this finite range. 

It is easy to verify that asymptotic values of the function B(x) at x — > — oo are directly related 
to those of the potential V which in this case plays the role of an effective cosmological constant: 



so that negative B(— oo) correspond to a de Sitter (dS) asymptotic, with B{— oo) = it is flat 
and with B{— oo) > it is anti-de Sitter (AdS). The solutions obtained may be classified by 
this asymptotic behavior and by the number and nature of horizons appearing there. The latter 
correspond to regular zeros of the function B(x) . It turns out that inclusion of the electromagnetic 
field makes the solutions much more diverse than it was found previously for purely scalar-vacuum 
configurations [1, 2, 17, 18] . 

3.2 Symmetric configurations 

To begin with, from (16) it follows that B(x) is an even function if and only if p = 0, hence m = 
(2/3)7rg 2 . Then V(x) is also an even function. Such symmetric configurations are asymptotically 
flat at both ends, x — > ±oo, and can be classified as follows (see the corresponding curves in Fig. 1): 

Al, A2: Twice asymptotically flat (M-M) wormholes. The curve A2 contains a minimum of B(x) 




(21) 



at x = . 



A3: Extremal regular black holes (M-M), with a double horizon (curve A2). 
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A4: Non-extremal regular black holes (M-M), with two simple horizons (curve A3). 

The abbreviation (M-M) stands here for two flat (Minkowski) asymptotic regions; we will also use 
similar notations for de Sitter (dS) and anti-de Sitter (AdS) asymptotic behaviors. 

The symmetric models form a one-parameter family, depending on q ; clearly, at q smaller than 
those appearing in Fig. 1 we also obtain wormholes (the simplest of them is with q = m = and 
V = 0, it is the Ellis massless wormhole [3,4]), while at larger q there are non-extremal regular 
black holes. The critical value of q that separates them is q ~ 0.825 (and m ~ 0.713), at which 
there emerges a double horizon corresponding to a double root of B(x), hence a regular extremal 
BH. 

It is of interest that in the narrow range of q in which the behavior of B(x) drastically changes, 
the potential V(x) changes very little. We also notice that at large \x\ the potential takes small 
positive values. It is not by chance since in the general case (20) V(x) behaves at large x as follows: 

Kw= §-i? +o(i ~ 7) - (22) 

3.3 Asymmetric configurations 

Concerning asymmetric configurations, it is natural to expect a critical behavior, i.e., transitions 
between different types of models, at values of m and q close to those appearing in Fig. 2 (but 
certainly with p ^ 0). This idea is confirmed by a direct inspection, and Fig. 2 (left) shows the 
corresponding five modes of the behavior of B(x) at m — 0.725: 

Bl: A black universe (M-dS) with a single simple horizon. 

B2: A black universe (M-dS) with two horizons (simple and double). 

B3: A black universe (M-dS) with three simple horizons. 

B4: A regular extremal black hole (M-AdS) with a double horizon, asymptotically AdS at the far 
end (x — > — oo). 

B5: A wormhole (M-AdS), asymptotically AdS at the far end (x — > — oo). 

The shape of the potential V(x) (Fig. 3, right) corresponds to Eq. (21): it is certainly zero at 
the flat asymptotic and is of the opposite sign to that of B at the other end. 

A somewhat different picture is observed if we slightly move down the mass and charge values, 
see Fig. 3 corresponding to m — 0.7. A qualitatively new feature as compared to Fig. 2 is that 
the function B(x) corresponding to a double horizon between two R- regions (curve C4) has a 
negative limit as x — > — oo. As a result, it is a black universe model instead of an M-AdS regular 
BH. Accordingly, the global causal structure characterized by the Carter-Penrose diagram is quite 
different (Fig. 3, bottom). 

Less diverse is the solution behavior at larger values of the parameters, as exemplified in Fig. 5: 

Dl, D2: Black universes (M-dS) with a single simple horizon, 
D3: Regular black holes (M-AdS) with two simple horizons. 
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Figure 2: Plots of B(x) (top, left) and V{x) (top, right) for asymmetric configurations at parameter 
values close to critical ones. The parameters are: m = 0.725 and q = 0.821, 0.8242, 0.829, 0.835, 0.839 
for curves B1-B5, respectively. In the corresponding Carter-Penrose diagrams, the letters R and T label 
R- and T-regions, respectively, while indices near R and T enumerate ranges with a certain sign of B(x) 
along the x axis from left to right. The symbols x n enumerate horizons (zeros of B(x)) from left to right: 

x\ < X2 < The diagram for B3 occupies the whole plane except the dashed triangles. The diagram for 

B2 cannot be placed in a plane since there are overlapping regions, as shown in the picture: the R3 region 
drawn below has a flat infinity (x — > 00) on the right, while in another R3 region, drawn turned up, such 
an infinity is on the left. 



The curve D3 corresponds to one more type of global causal structure: the Carter-Penrose 
diagram (Fig. 4, right) is the same as for a non-extremal Reissner- Nordstrom BH, but instead of a 
Reissner-Nordstrom central singularity we have an AdS infinity. 

So far we have been assuming that the space-time is asymptotically flat as x — > 00. It is clear 
that if we abandon this assumption, then the number of possible qualitatively different globally 
regular configurations in the scalar-electrovacuum system under consideration will be still larger. 
To see how they can look, let us note that in Eq. (16) the constant B is additive. Therefore, 
changing B , we simply move up or down the plot of B(x), thus changing the asymptotic behavior 
and the number and nature of horizons in our model. For instance, if we slightly move down the 
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Figure 3: Plots of B(x) (top, left) and V(x) (top, right) for asymmetric configurations with the parameter 
values m = 0.7 and q = 0.806, 0.8086, 0.811, 0.8148 for curves C1-C4, respectively. In the Carter-Penrose 
diagram for C4 shown below, the notations are the same as in Fig. 2. 




Figure 4: Plots of B(x) (left) and V(x) (middle) for asymmetric configurations with the parameter values 
m = 1 and q = 0.968 (curve Dl), q = 0.974 (curve D2) and q = 0.98 (curve D3). The right panel shows 
the Carter-Penrose diagram for D3. 

curve A3 in Fig. 2, we obtain a configuration with two de Sitter asymptotics (dS-dS), separated by 
four simple horizons, see Fig. 5. 

In the same way it is easy to obtain a number of other configurations with dS and AdS asymptotic 
behaviors. 
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Figure 5: A symmetric configuration (dS-dS) with four simple horizons, m = 0.739, q = 0.84, Bq = 
— 37rm/2 + 7r 2 (7 2 /4 — 0.05. Plots of B{x) (top, left) and V(x) (top, right) and the Carter-Penrose diagram 
(bottom), with the same notations as before. The diagram occupies the whole plane except the segments 
shown by wavy lines (de Sitter infinities). 

4 Discussion 

Scalar-vacuum configurations with a self-interacting phantom scalar field have been considered 
in [1,2] (see also references therein); they included M-M and M-AdS wormholes and black universes. 
In the present paper, we have obtained similar models with an electromagnetic field added and 
found that its inclusion leads to a greater diversity of qualitatively different configurations. More 
specifically, even being restricted to solutions which are asymptotically flat as x — > oo and have 
m > 0, we have found as many as 10 types of models, classified by the types of asymptotic behavior 
and the number and nature of horizons. At zero charge q we return to the situation discussed in [1,2], 
with only two configuration types: M-M wormholes with m = (its analogue is represented here 
by the curve Al) and black universes with a single simple horizon (a similar behavior of B(x) is 
shown here, e.g., by the curves Bl and CI); also, M-AdS wormholes were obtained there but only for 
m < . As already mentioned, the reason for such a narrow choice is that in a pure scalar- vacuum 
system the field equations forbid the function B(u) to have a regular minimum [16]. 

Different types of regular configurations obtained here, which are asymptotically flat as x — > oo 
and have a nonnegative Schwarzschild mass, are summarized in the table. 

The general landscape of our solutions with a flat infinity on the right end and a positive 
Schwarzschild mass can be characterized by a map drawn in the (m,q) plane (Fig. 6). It can be 
seen that the most generic are M-AdS wormhole solutions and black universes with a single simple 
horizon: they exist for all values of m and are actually the same types of solutions that have been 
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Table 1: Types of asymptotically flat solutions with m > 



Solution type 
(curve number) 


Configuration type, asymptotics 
(x — > +oo) — (x — > — oo) 


Horizons: number, order n 
[disposition of R- and T-regions] 


A1,A2 


M - M wormhole 


none [R] 


A3 


M - M extremal black hole 


1 hor., n = 2 [RR] 


A4 


M - M black hole 


2 hor., n = 1 (both) [RTR] 


Bl, CI, C4, Dl, D2 


M - dS black universe 


1 hor., n = 1 [TR] 


B2, C2 


M - dS black universe 


2 hor., n = 2 and n = 1 [TTR] 


C4 


M - dS black universe 


2 hor., n = 1 and n = 2 [TRR] 


B3, C3 


M - dS black universe 


3 hor., n = 1 (for all), [TRTR] 


D3 


M - AdS black hole 


2 hor., n = 1 [RTR] 


B4 


M - AdS extremal black hole 


1 hor., n = 2 [RR] 


B5 


M - AdS wormhole 


none [R] 




m 

Figure 6: Map of asymptotically flat solutions with m > 0. Right panel: a domain of interest enlarged 

obtained before with q = [1]. One more generic type of models is formed by regular BHs with 
two horizons, which appear only with sufficiently large charges. Solutions with extremal horizons 
appear on separatrices between the main domains on the plane, while solutions with three simple 
horizons are also generic since they occupy a certain area on the plane, but this area is actually a 
very narrow band, it is almost invisible if we do not specially adjust the scale, see the right panel 
in Fig. 6. 

We conclude that the present field system creates quite a number of diverse models, making us 
substantially widen the list of possible regular BH configurations as compared, e.g., with [2]. Such 
models can be of interest both as descriptions of local objects (black holes, wormholes) and as a 
basis for building singularity-free cosmological scenarios. An important feature of such cosmologies, 
different from the great majority of nonsingular models described in the literature, is that the 
cosmological expansion starts from a Killing horizon (this phenomenon can be termed a Null Big 
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Bang [19,20]) beyond which, in the absolute past, there is an asymptotically flat static region. There 
is another kind of configurations with a Null Big Bang where a static region, instead, contains a 
regular center [19,21]; as in the present paper, the models described there can possess multiple 
horizons and have a de Sitter asymptotic behavior at late times. 

An important point concerns the value of the global magnetic field that exists in our Universe if 
it can be described by a model more or less like ours. Let us use the (probably) most conservative 
estimate, according to which a lower limit on the magnetic field strength is B > 1CT 18 G [22]. On the 
other hand, the present scale factor a ~ 10 28 cm approximately coincides with the quantity r{u) in 
the metric (2). Therefore, since B oc r~ 2 [see (6)], we can roughly estimate the field values at earlier 
stages of the evolution. For instance, at recombination (a/ do ~ 10 -3 ), when the electromagnetic 
radiation decoupled from matter, the magnetic field was still weak enough, B ~ 10~ 12 G, but at 
the stage of baryogenesis (a/a ~ 1CT 12 ) it was of the order of 10 6 G. 

Observations of the cosmic microwave background (CNB) show that our Universe is highly 
isotropic: at recombination, the degree of anisotropy did not exceed 10~ 6 . If this condition holds, 
spherically symmetric models like ours (belonging to the Kantowski-Sachs class), being anisotropic 
by construction, can still conform to observations [23]. 

The above condition constrains the global magnetic field strength allowed by the observed CMB 
isotropy. The CMB energy density pcmb and that of the the magnetic field, p ma gn, are both 
proportional to r -4 as long as the Universe is approximately isotropic, hence their ratio is constant 
and is the same at recombination and at present. But at present Pcmb ~ 10~ 33 g/cm 3 , hence we 
should require p magn < 1CT 39 g/cm 3 , which in turn means \B\ < 1CT 8 G. We see that this condition 
is easily satisfied by the fields under consideration. 

An upper limit on the classical magnetic field description seems to follow from the work of Am- 
bjorn and Olesen [24] who pointed out that the Weinberg-Salam model of electroweak interactions 
shows an instability at B > 10 24 G, connected with emergence of a tachyonic mode. If at present 
\B\ ~ 10~ 18 G, then the maximum admissible value of B corresponds to r ~ 10 7 cm = 100 km - 
it is the minimum admissible value of b = minr(w) in our models. 

The possible viability of models like those considered in this paper depends on their stability 
under various kinds of perturbations. Most of the known scalar-vacuum wormhole and black- 
universe solutions proved to be unstable under radial perturbations [25-27], and it is of interest to 
find out whether or not they can be stabilized by electric or magnetic fields. We hope to consider 
this problem in the near future. 
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